Nonlocal memory effects in the dynamics of open quantum systems 
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We explore the possibility to generate nonlocal dynamical maps of an open quantum system 
through local system-environment interactions. Employing a generic decoherence process induced 
by a local interaction Hamiltonian, we show that initial correlations in a composite environment can 
lead to a nonlocal open system dynamics which exhibits strong memory effects although the local 
dynamics is Markovian. In a model of two entangled photons interacting with two dephasing envi- 
ronments we find a direct connection between the degree of memory effects and the amount of cor- 
relations in the initial environmental state. The results demonstrate that, contrary to conventional 
wisdom, enlarging an open system can change the dynamics from Markovian to non-Markovian. 
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Coupling a quantum mechanical systems to an external 
environment causes the system to lose information to its 
surroundings. Since merely all realistic quantum systems 
are open, understanding and controlling the dynamics 
arising from the presence of the environment is of central 
importance in present-day research [TJ [2] . The standard 
approach to the dynamics of open quantum systems em- 
ploys the concept of a quantum Markov process which 
is given by a semigroup of completely positive dynam- 
ical maps and a corresponding quantum master equa- 
tion with a generator in Lindblad form [21 H]- However, 
many quantum systems exhibit non-Markovian behavior 
in which there is a flow of information from the environ- 
ment back to the open system, signifying the presence of 
quantum memory effects [StilOj. 

For many processes occurring in nature the approxi- 
mations allowing a simple Markovian description are not 
applicable. It is known, for example, that strong system- 
environment couplings, structured and finite reservoirs, 
low temperatures, as well as the presence of large initial 
system-environment correlations can give rise to mem- 
ory effects in the open system dynamics. The recogni- 
tion of the importance of non-Markovian processes has 
initiated many essential steps towards the development 
of a general consistent theory of non-Markovian quantum 
dynamics [TllfTT] as well as achievements in the experi- 
mental detection and control of memory effects [HI US] . 



In this Letter we introduce a hitherto unexplored 
source for quantum memory effects, namely the pres- 
ence of initial correlations between the subsystems of a 
composite environment which interact locally with the 
subsystems of a composite open system. It is demon- 
strated that correlations between the environmental sub- 
systems can generate a nonlocal quantum process from 
a perfectly local interaction Hamiltonian. We will show 
further that a nonlocal decoherence process can lead to 
non-Markovian behavior although the local dynamics of 



both subsystems is Markovian. These features are dis- 
cussed by employing two theoretical models, namely a 
generic decoherence model of two qubits interacting with 
correlated multimode fields, and an experimentally real- 
izable model of entangled down converted photons trav- 
eling through birefringent media. We thus find a new, 
experimentally controllable source for memory effects in 
a quantum dynamical process. Besides the practical im- 
portance of the result in the physical realization and con- 
trol of dynamical processes, it also reveals an unexpected 
feature about the nature of non-Markovian dynamics of 
composite quantum systems: Enlarging the open system 
can actually turn the dynamics from a Markovian to a 
non-Markovian regime. 

We consider an open system S consisting of two sub- 
systems labeled by an index i = 1,2, and an environ- 
ment E which is also composed of two subsystems. We 
assume that there are only local system-environment in- 
teractions, i.e., that subsystem i of 5 interacts only with 
its environment i of E. The local interactions are de- 
scribed by unitaries Ui(t), and S and E are supposed 
to be uncorrelated at the initial time. The open system 
state at time t is thus given by 



V(o)) 



where $12 (t) represents the quantum dynamical map de- 
scribing the time evolution of S. If the two environments 
are initially uncorrelated, p^Ei^) — Pe^) ® t^is 
map factorizes and the dynamics of S is given by a prod- 
uct of local maps, ^nit) = 4>i(t)(g)$2(i)- However, when 
Pil(O) exhibits correlations $12 (i) does not in general fac- 
torize and the environmental correlations may give rise to 
a nonlocal process even when the interaction Hamiltonian 
is purely local. For a local map the dynamical properties 
of the subsystems completely determine the global sys- 
tem dynamics, but when the map is nonlocal the global 
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system can exhibit features which are not present in the 
dynamics of the individual subsystems. Here we explore 
especially quantum memory effects arising from a nonlo- 
cal dynamics. 

We consider a dephasing map for two qubits of the 
general form 



ba*K'^{t) 
ca*Kl{t) 



ab* K2{t) ac*Ki(t) ad* Ki2{t) \ 
bc*Ai2{t) bd*Ki(t) 



cb*Kl2{t) \cY 



\da*i<{2{i) db*K\{t) dc*K^(i) 



cd*K2{t) 

\d\' 



I 

(1) 

where the initial state of the two qubit system is a pure 
state given by 



IV'12) =a|00)+6|01>+c|10>+d|ll) 



(2) 



The corresponding dynamics for subsystems 1 and 2 are 
given by p\;{t) = tr2 [pf{t)\ and p|(t) = tn [pf{t)\. 
The states p\{t) and are fully determined by the 

functions ni{t) and ^2(^)1 depending neither on Ki2(t) 
nor on Ki2{t). The interaction Hamiltonian is assumed 
to be local, i.e., we have 



Hint{t)^Xlit)Hi+X2{t)H2, 



(3) 



where the function Xiit) is 1 for < t < and 
zero otherwise. Here, tf and tl denote the times the 
interaction is switched on and switched off in system 
i, respectively. Since the local Hamiltonians Hi com- 
mute, the time evolution of the total system is given by 



exp 



l*(0)). 



We will further 
t 



-I j;dt'H,^,{t' 

denote the local interaction times as ti(t) — Xi{t')dt' 
and for convenience we will not explicitly write the time 
dependence of ti. 

Before turning to the details of the physical systems 
under study, let us briefly discuss the concept of mem- 
ory effects. Memory effects are quantified in T3 by em- 
ploying the trace distance D{pA, Pb) = \^Apa — Pb\ be- 
tween two quantum states pA and pB- This quantity 
can be interpreted as a measure for the distinguishability 
of the two states [20 22^ . In view of this interpretation 
the characteristic feature of a non-Markovian quantum 
process is the increase of the distinguishability, i.e. a re- 
versed flow of information from the environment back 
to the open system. Through this recycling of informa- 
tion the earlier states of the open system influence its 
later states [S] , which expresses the emergence of mem- 
ory effects in the open system dynamics. The measure 
for non-Markovianity is written as 



A/'($) = max 

Pa.b(O) 



dta{t,pA^B{0)), 



(4) 



CT>0 



where a{t, pa,b{0)) — jiD{pA{t), PB{t))- Here, the time 
integration is extended over all subintervals of time in 
which the rate of change of the trace distance a is pos- 
itive, and the maximum is taken over all pairs of initial 



states. The quantity in Eq. Q thus measures the maxi- 
mal total amount of information which flows from the en- 
vironment back to the open system over the whole time 
evolution. 

First, we study a generic model of two qubits inter- 
acting with correlated multimode fields. The local inter- 
action Hamiltonians of Eq. ^ are Hi — J2k '^lidkb^k + 
9k^k)- We assume that the interaction strengths in both 
systems are identical, gl = g^- The local time evolution 
of the systems is then given by the unitary 

U.{t) = eKp{alY,{b^^k{t,) - blCkik))}, (5) 
fc 

where £,k{ti) — gfc(l — e"^''*')/wfc. The local unitary of 
Eq. ^ acts in the following way: 



U,{t)\0)®\rj) = |0) 



k 

where -D(^fe) is the displacement operator for the fcth 
mode. Let us take as the initial state |^'(0)) = |'0i2) ® 
1 7712), where IV'12) is given by Eq. ^ and |?7i2) = 
®fc l^l?)- "^^"^ decoherence process is then given by 
Eq. (0, where = {vl°2KA), ^2{t) = {rillK^), 

«i2(t) = {vlWrHD, and Ai2(t) = {vllKD with 



1^12" W) = (g) [D{{~ir+'Uh)) ® D{{-ir+^Ut2))] 1^712) 
fc 



After some algebra, one finds 
Ki(i) 

K2(i) 

^12(0 

A12W 



nXfc(-2efe(ti),0), 

k 

nXfc(0,-2a(i2)), 

k 

nxfc(-2efc(ti),-2efc(t2)), 

k 

J|xfc(-2efc(ti),+2a(t2)), 



where Xk{x,y) is the characteristic function of \vii2)- 

Let us consider a two-mode Gaussian state with the 
characteristic function Xk{x,y) = Xfc(Ai, A2, A3, A4) = 
exp(— iA-^crA), where Ai — 5R[x], A2 = A3 = '^[y]^ 

A4 = and 



A C 
B 



is the covariance matrix of the state. Let us take A = 
B — I and C = cl. Now the state is uncorrelated iff 
c = 0. We can write 



Xkix.y) = exp 



+ \y\^ + c{xy* +x*y)} 
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For c = — 1 we get Xk{x,y) = exp [— ^ja; — yp] and 
Ki2(i) = exp [-2J2k l^kih) - 6(^2)1^]- Performing the 
continuum limit we obtain 



Ki2(t) = exp<^ -4 



where Jiuj) is the spectral density of the reservoir. For 
an ohmic spectral density J(w) — aco exp(—uj/ujc) with 
coupling constant a and frequency cutoff lUc we have 



A12W 



'Ait)) 



,2.2 



m) 



-2a 



-2a 



{l + col\h{t)-Ut)\') 

Kl{t)Kl{t)/Ki2{t). 



The maximization over the initial states in Eq. Q for 
different values of c is presented in Fig. (TJa). The trace 
distance dynamics of the subsystems 1 and 2 as well as 
the global dynamics are presented in Fig. [2] We clearly 
see that the trace distance in the subsystems 1 and 2 con- 
tinuously decreases, but for the total system the trace dis- 
tance does indeed increase: We obtain a dynamics which 
is locally Markovian but globally exhibits nonlocal mem- 
ory effects. 




FIG. 1: (color online) Maximization over the pairs of initial 
states. The (blue) dots represent the increase of the trace 
distance for 1000 randomly drawn initial pairs of states, (a) 
Two qubits interacting with multimode fields. The measure 
for non-Markovianity for different values of c and for a — 1. 
The (red) triangles represent the measure for the maximizing 
pair ^(|00)±lll)), and the (red) squares for ^(|01) ± 1 10)). 
(b) Two photons moving through quartz plates. The measure 
for non-Markovianity for different values of the correlation 
coefficient K, and a fixed cK^TAn = 1. The (red) triangles 
represent the measure for the maximizing pair -^(|_ff//) ± 

\VV}), and the (red) squares for -^{\HV} ± IVH}). 

As our second example we examine an experimentally 
realizable model of a pair of entangled photons subjected 
to local birefringent environments [531 HI]- The pho- 
ton pair is created in a spontaneous parametric down- 
conversion process after which the photons separate, 
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FIG. 2: (color online) The trace distance dynamics for the two 
qubits interacting with correlated multimode fields. We take 
a = 1, = 0, = 1 = t| and = 2 in units of w<r^ The 
solid (blue) lines represent the trace distance with different 
values of c for the global dynamics of the two qubits for the 
maximizing pair of initial states. The dashed (red) line and 
the dotted (green) line give the trace distance evolution for the 
initial states l/\/2(|0) ± |1)) in systems 1 and 2, respectively. 



traveling along different arms i = 1,2, and move through 
different quartz plates. When a photon enters a quartz 
plate a local interaction between the polarization degrees 
of freedom (forming the open system) and the frequency 
degrees of freedom (forming the environment) is switched 
on. The Hamiltonian describing the local interaction in 
Eq. ([3| of a photon induced by the corresponding quartz 
plate is given by 



nv\V){V\+nH\H){H\ <^\uj,){u, 



where |A) ® denotes the state of a photon in arm i 
with polarization \ = H^V (horizontal or vertical) and 
frequency w^. The refraction index of the polarization 
state A is denoted by nx. The total initial state is given 
by 1^(0)) = IV'12) ® J duJiduj2 g{^i,0J2) \'^i,t^2) , where 
1^12) = a \HH) + b \HV) + c \VH) + d \VV). Initially 
the environment formed by the mode degrees of freedom 
is thus in a correlated state with g{uji,Lij2) denoting the 
amplitude of finding a photon with frequency loi in arm 
1 and a photon with frequency uj2 in arm 2. The cor- 
responding joint probability distribution will be denoted 
by P{uji,uj2) = |g(wi,a;2)p. 

The state of the open system (polarization states) 
at time t is of the form of Eq. ([T]) with the func- 
tions Ki(t) = G(Anti,0), K2{t) = G{0,Ant2), Ki2(t) = 
G{Anti,Ant2), and Ai2{t) = G{Anti,-Ant2), where 

G(ti,T2)= / dL^idco2P(wi,c^2)e-'("^"^+"="=^ 



is the Fourier transform of P(a;i, ^2) and An = ny — nn 
is the birefringence. Note that although the Hamilto- 
nian of (|3| is a sum of local interaction terms, the cor- 
responding dynamical map $12(0 is ^ product of local 
dynamical maps if and only if Ki2(t) = Ki(t)K2it) and 
Ai2(i) = '«i(t)K2(t). This is the case only when the joint 
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frequency distribution P(uji,uj2) factorizes, i.e., when the 
frequencies uii and W2 are uncorrelated. 

To characterize the correlations in the initial environ- 
mental state we introduce the covariance matrix C = 
{Cij) with elements Cij = {tOiUJj) — {uji){ujj). We will 
assume that both the means and the variances of loi and 
UJ2 are equal, i.e., (cji) = {L02) — ujf)l2 and Cn = C22 — 



To quantify the frequency correlations 



we use the correlation coefficient K = C12/ ^/CuC^^ ~ 
C12/C11. We have \K\ < 1, where the equality sign holds, 
i.e., K = ±1, if and only if wi and UJ2 are linearly related. 
Let us take a Gaussian frequency distribution 



P{uJi,UJ2) 



1 



27rVdetC 



-i(c3-M)^C-i(cD-(<:)) 



(6) 



where cD ~ {uji,uj2)^ and (u) = {{oji), {uj2))'^ ■ One can 
easily find the Fourier transform of this distribution. 

We assume for simplicity that the total interaction 
times for both photons are equal, denoting it by T = 
t|, and that the quartz plates are mounted 



4-S _ 4-t 
Li In 



one after another, i.e., t\ — i|. We can then derive 
an analytic expression for the measure ^ of the non- 
Markovianity of the process. The maximization over the 
pair of initial states in Eq. Q is illustrated in Fig. [TJb). 
Using Eqs. ([T]) and ([t]) we obtain the time dependence 
of the trace distance for the maximizing initial pairs. 



D{t) = exp 



Cu{tl+tl-2\K\tit2) 



During the interaction of the photon in arm 1 the trace 



distance first decreases from 
value Di = exp 



he initial value 1 to the 
2 -CiiT"^ . The subsequent inter- 
action of the photon in arm 2 depends on the function 
f{t2) = 4- 2\K\Tt2 for t2 e [0,r]. This function de- 
creases monotonically in the interval [0, \K\T] from the 
value /(O) = to the value f{\K\T) = -{KTf, which 
means that the trace distance increases over this inter- 



val to the value D2 = exp 



An' 



It 



follows that the non-Markovianity measure is given by 

M - D2-Di = e-Kii(A"T)^ UCi.iAnTfK- _ i] . (g) 

This equation establishes a direct connection between the 
measure for non-Markovianity and the degree of corre- 
lations in the initial environmental state as quantified 
by the correlation coefficient K. We also see that the 
process is Markovian if and only if if = 0. The rela- 
tion ([s]) is further illustrated in Fig. [sj where we have 
plotted the frequency distribution P(a;i,a;2) and the dy- 
namics of the trace distance for three different values of 
the correlation coefficient, K = 0.0, —0.5, —1.0 (anticor- 
relation). One clearly observes that when the frequen- 
cies wi and UJ2 become more anticorrelated the dynamics 




FIG. 3: (color online) The connection between the frequency 
distribution P(a;i,cJ2) (top) and the dynamics of the trace 
distance (bottom) for different values of K. (a) K = —1, (b) 
K = -0.5, (c) K = 0. The unit of time is y/U^TAn. 



becomes more non-Markovian. In general, we conclude 
that the reduced dynamics of the two-photon polariza- 
tion state is non-Markovian whenever the frequency dis- 
tribution P[uji,uj2) exhibits correlations. This behavior 
occurs globally, i.e., when we study the dynamics of the 
composite state of both photons. However, if one ob- 
serves the local dynamics of either of the photons the 
process is always Markovian. 

Summarizing, we have found a new source for mem- 
ory effects in the dynamics of open quantum systems. 
We studied a generic dephasing model as well as a the- 
oretical scheme which is experimentally realizable with 
current technology. The initial correlations within the 
environment induce a nonlocal map which gives rise to 
memory effects in the open system dynamics. Locally, 
each subsystem follows Markovian dynamics, but glob- 
ally they are exposed to memory effects even though 
the interaction Hamiltonian is local. Since for classical 
stochastic processes a non-Markovian process can be em- 
bedded in a Markovian one by a suitable enlargement 
of the number of relevant variables, the general view 
has been that enlarging the open quantum system under 
study tunes the dynamics towards a Markovian behav- 
ior. This can be done for certain non-Markovian quan- 
tum processes as well j25U28j . but we see that also the 
exactly opposite behavior can occur, i.e., enlarging the 
subsystem can bring the dynamics from a Markovian to 
a non-Markovian regime. 
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